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We study the structures of the unstable Be isotopes with the theoretical method of 
CN| ! antisymmetrized molecular dynamics. It is found that various structures of the excited 
states appear in the low-energy region of neutron-rich Be nuclei. Focusing on the 2a 
clustering, we analyze the intrinsic structures with the help of the experimental data of 
Gamow- Teller transitions. 

b 

a ■ 1. INTRODUCTION 



Recently, the study of unstable nuclei is one of the hot subjects. In the light unstable 
nuclei, the experimental informations near the drip lines increase rapidly and many inter- 
esting phenomena are found. The exotic structures such as the halo, the skin structures, 
and the clustering structures are suggested in some light unstable nuclei. Are the features 
of the clusters in unstable nuclei similar to those in the stable nuclei ? The shift of the 
magic number and the deformations of unstable nuclei are the concerned problems. 

Our aim is to research on the excited states of unstable and to study systematically the 
various structures of unstable nuclei. We apply a theoretical method of antisymmetrized 
molecular dynamics (AMD). The method of variation after projection(VAP) in the frame- 
work of AMD has been already proved to be useful to investigate the excited states of 
light unstable nuclei |],|2| as well as stable nuclei 0. AMD is a powerful approach which 
is applicable to the excited states of general nuclei. The wave function of AMD can rep- 
resent various clustering structures as well as shell-model-like structures, because there is 
no inert core in the framework. 

In the previous study on the excited states of 10 Be with AMD MM, we found that 
the molecule-like structures may coexist in the low-energy region with the ordinary shell- 
model-like states. The data of the transitions such as E2 and Gamow- Teller transitions 
are so helpful to classify the excited levels into rotational bands. The reader is referred 
to the paper M for the detailed discussions. 
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2. FORMULATION 

In this section we explain the formulation of AMD for the study of the nuclear structure 
of the excited states. 

2.1. Wave function 

The wave function of a system is written by AMD wave functions, 

® = c& amd + (!& AM d + -'-- (!) 
An AMD wavefunction of a nucleus with mass number A is a Slater determinant of 
Gaussian wave packets; 

$AMD&) = -jt.A{<Px,<P2,---,<PA}, (2) 



<Pi = foiXfcTj : < 



(3) 



( 1-1 ; ' 

where the zth single-particle wave function ^ is a product of the spatial wave function, 
the intrinsic spin function and the iso-spin function. The spatial part is presented by 
complex parameters Xu, X21, X&, of the centers of Gaussians, and x& is the intrinsic spin 
function parameterized by £j. T{ is the iso-spin function which is fixed to be up (proton) or 
down(neutron) in the present calculations. Thus an AMD wave function is parameterized 
by a set of complex parameters Z = {X ni ,^i} (n — 1,3 and i — 1, A) which determine 
the centers of Gaussians of the spatial parts and the directions of the intrinsic spins for 
all single-particle wave functions. 

If we consider a spin-parity eigen state projected from a AMD wave function, the total 
wave function is linear combinations of Slater determinants, 

$(Z) = Pfiic>*AMD(Z) = (1 ± P) J dQD{; K ,(Q)R(Q)<l> AMD (Z), (4) 

where P is a parity projection operator. The integrations for the Euler angles are numer- 
ically calculated by a summation of mesh points on Q. 

In principal the total wave function can be a superposition of independent AMD wave 
functions, 

$ = cP^amdCZ) + c'P^AMniZ') + ■■■■ (5) 
2.2. Energy variation 

We make variational calculations to find the state which minimizes the energy of the 
system; ($|//|$)/($| ( I ) ) by the method of frictional cooling one of the imaginary 

time methods. The time development of the wave function $(Z) is given by the equations, 
dX nk , 1 d ($(Z)|ff|$(Z)) 

-dT = {x + %fi) ihdxT k <$(z)|$(z)> ' (n = 1 ' 3 fc = 1 ' A) (6) 
de* a , ■ x 1 ^ ($(Z)|j/|$(Z)) 

^ = (A + z/i) ^ (^(z)|<&(z)) ' (A; = 1 ' A) (7) 

with arbitrary real numbers A and /i < 0. It is easily proved that the energy of the 
system decreases with time. After sufficient time steps of cooling, we obtain the optimum 
parameters for the minimum-energy state. 
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2.3. Lowest J^ 1 states 

In order to obtain the wave function for the lowest state, we perform the energy 
variation for the spin parity eigenstates projected from an AMD wave function. In this 
case the trial function is $ = P'l^ K ,^ amd^) ■ For the preparation of VAP calculations, 
first we make variational calculations before the spin projection in order to choose an 
initial wave function for each parity and an appropriate K' quantum for each spin-parity 
state. Then we perform VAP calculation for 

(P^±,^ MO (Z)|i?|P J ^,$^(Z)}/(P^(Z)|P^±,(Z)) (8) 

with the adopted K' quantum from the initial state. In the VAP procedure, the 3-axis 
of Euler angle in the total-angular momentum projection is not necessarily same as the 
principal z-axis. The direction of the intrinsic deformation in the body-fixed frame is 
automatically determined in the energy variation depending on the adopted K' quantum. 
However, in many cases, the approximately principal axis z obtained by a VAP calculation 
is found to be almost equal to the 3-axis. 

2.4. Higher excited states 

As mentioned above, with the VAP calculation for $(Z) = P^^amdC^), we obtain 
the set of parameters Z = Z^ which presents the wave function for the first J ± state. To 
search the parameters Z for the higher excited J states, we make the energy variation 
for the orthogonal component of a spin-parity projected AMD wave function to the lower 
states by superposing the wave functions. If the AMD wave functions for the higher 
excited states have energetically local minimums, we do not need orthogonalization in 
VAP calculations. 

2.5. Expectation values 

After VAP calculations for various states, we obtain many intrinsic states <&\md-> 
^amdi' ' 'i ^™mo' Finally we improve the wave functions by diagonalizing the Hamilto- 
nian matrix {P J Jk'^AMD\ H \ P Mk"^AMD) and the norm matrix (Pmk'$WdI p m^» $ aa/d) 
simultaneously with regard to for all the intrinsic states and (K',K"). In compari- 
son with the experimental data, the theoretical values are calculated with the final states 
after the diagonalization. 

3. RESULTS 

We calculate the excited states of 12 Be and 11 Be. The adopted interactions are the MV1 
case3 (the finite two-body and the zero-range three-body interactions for the central 
force, and G3RS for the spin-orbit force || . The coulomb interaction is approximated by 
a sum of seven Gaussians. In the present calculations, we adopt the Majorana parameter 
m = 0.65 and omit the Bartlett and Heisenberg terms in the two-body central force. As 
for the strength of the spin-orbit force, we try two cases of the parameter, m; s =3700 MeV 
and ui s =2500 MeV. 

The resonance states are treated as a bound-state approximation in AMD, because 
the restriction of the Gaussian forms of the wave functions makes an artificial barrier for 
outgoing particles. 
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Figure 1. Energy levels and transitions of 12 Be. The adopted parameters for the calcula- 
tions are Majorana parameter m = 0.65 and the strength of the spin-orbit force u\ s = 3700 
MeV. 
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3.1. 12 Be 

The energy levels of 12 Be are displayed in Figure |T]. The theoretical results are calculated 
with the stronger spin-orbit force as Ui s = 3700 MeV which can reproduce the abnormal 
parity of the ground l/2 + state in 11 Be. Many excited states are seen in the low-energy 
region. Even though the 12 Be has the neutron magic number 8, it is surprising that 
the calculated ground state is not the ordinary state with the neutron closed p-shell, 
but a 2 particle 2 hole state with the developed clustering structure. As the evidence 
of the breaking of the closed p-shell, the calculations reproduce well the small value of 
the Gamow-Teller(GT) transition strength from 12 Be(0 + ) to 12 B(1 + ). The experimental 
B(GT) is about one third as small as the value expected for the closed p-shell 12 Be(0 + ). 
In other words, the weak GT transition is one of the proofs for the disappearance of the 
neutron magic number 8 in neutron-rich Be isotopes. It is consistent with the shell model 
analysis by T. Suzuki et al. ||. We also find many excited states with the developed 
molecule-like structures which construct the rotational bands K n = Of, = K n = 
Of in the low-energy region. Since the ground state has the large deformation in the 
present calculations, The predicted B(E2; 2f — > 0^) in the ground band is 14 e 2 fm 4 and 
is much larger than the value 7 e 2 fm 4 of B[E2; 2f — > Of) in the second K n = + band 
which consists of the closed p-shell states. 
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Figure 2. Energy levels of n Be. The adopted parameter for the strength of the spin-orbit 
force is ui s = 3700 MeV for case(l) and ui s = 2500 MeV for case(2). 
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3.2. n Be 

In this subsection, we discuss the excited states of n Be. In Fig 0, the calculated energy 
levels are shown with the experimental data. We use the rather strong spin-orbit force in 
case (1) and the weaker LS force in case (2). The binding energies are underestimated as 
58 MeV in case (1) and as 54 MeV in case (2). It is not difficult to fit to the experimental 
binding energy by choosing a smaller value for the Majorana parameter. 

In the case (1) with the stronger LS force, the abnormal spin-parity l/2 + state comes 
down lower than the normal negative parity states, which is consistent with the experi- 
mental data. In both cases, we found some rotational bands K n = l/2~, K 71 = l/2 + and 
K n = 3/2~, although the order of the levels depends on the adopted interaction. The 
K w = 3/2~ band has a large moment of inertia and reaches to the high spin states. The 
highest spin is 15/2 - in case(l) and 13/2~ in case (2). It is expected that the K 71 = 3/2~ 
band may reach to the further high spin such as a 17/2 - state in the calculations with a 
smaller value of the Majorana parameter. 

Next we analyze the intrinsic wave functions of the obtained rotational bands. In the 
results of 11 Be, 2a clustering structures appear very often in many states. The schematic 
figures are presented in Fig. ||. Two alpha cores are found in the states of the lowest 
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Figure 3. Schematic figures for the intrinsic structures of 11 Be. (a), (b), (c), (d) correspond 
to K = 1/2", K = 1/2+, K = 3/2-, 3/2~(« 10 MeV), respectively. 






(a)K=l/2 _ (b) K=l/2 + (c) K=3/2~ (d)3/2" lOMeV) 



negative parity band K n = l/2~ (Fig||(a)). In the positive parity states, one neutron 
is in the sd orbit and the 2a clustering develops (Fig.^(b)). We think that one of the 
reason of the energy gain of the positive parity states is considered to be the clustering 
development. The interesting point is that the K n = 3/2~ band comes from the well- 
developed molecule-like structure with 2 neutrons in the sd orbit (Fig.|3](c)). The large 
moment of inertia results from the largely deformed intrinsic structure. At about 10 
MeV, we found another characteristic state where one of the 2 a-clusters completely 
breaks (Fig.|3|(d)). In this state the intrinsic spins of 2 protons couple up to be totally a 
unit. As mentioned above, various structures coexist in the low energy region in n Be. It 
is expected that the valence neutrons play important roles in these excited states. 

We discuss the Gamow-Teller(GT) transitions in order to analyze the level structures of 
n Be. First we consider /3 decay from n B. The experimental data of GT-type transitions 
from the ground state of n B are deduced from the 0° cross sections of the charge exchange 
reactions [|K|. Comparing the calculated log ft values with experimental data (Fig.|j), the 
observed three levels well correspond to the excited states in the lowest K n = 1/2" band. 
It is natural that the GT transitions from n B to these levels are reasonably strong because 
both the parent and the daughter states are ordinary p-shell states. 

The beta decays from n Li are shown in Fig.[|. We notice that the ideal 2a-cluster 
structures in the daughter n Be states forbid any GT transitions from n Li because of the 
Pauli forbidden. In other words, the GT transitions from 11 Li are allowed only if the 
daughter states of 11 Be have the component with the breaking of the a clusters. The 
calculations well agree with the experimental log ft values. According to the calculated 
results, the lowest 3/2~ and 5/2~ states in the K w = l/2~ band have the more than 10 
% breaking of the clusters. One of the reason for the weak GT transition to the lowest 
l/2~ state is the clustering structure of n Be. Another reason is suggested by Suzuki et.al 
[d as the halo structure of the parent 11 Li state, which may not be described enough 
in the present AMD calculations. In the calculations, transitions for the higher excited 
states in the K 71 = 3/2~ band are suppressed because of the a clusters in the developed 
molecule-like structures. The strong transition to the excited state at 8 MeV well fits to 
the calculated 3/2~ state (~10 MeV ) where one of the two a clusters completely breaks. 

As for the El transitions, the B(E1; l/2~ — > l/2 + ) of the present calculations is much 
smaller than the experimental data. We should describe the wave function more precisely 
to reproduce the strength B(E1). 
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Figure 4. Beta decays to the excited states of Be. GT transitions are calculated 
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4. SUMMARY 

We studied the structure of the excited states of 12 Be and n Be and analyzed the level 
structure with the help of the data of Gamow- Teller transitions. 

In 12 Be, it is found that the ground state has the developed clustering structure instead 
of the neutron p-shell closed state. The main component of the ground state is the 2p- 
2h state in terms of a simple shell model. That is to say the neutron magic number 8 
disappears in 12 Be. It is consistent with the small GT strength to 12 Be. Many molecule- 
like states make rotational bands in low energy regions. 

In the case of n Be, various structures are seen in the excited states. In the results 
rotational bands with the the developed clustering states are suggested. We discuss the 
Gamow- Teller transitions comparing the AMD calculations with the experimental data. 
The experimental data of GT transitions from 11 B well agree with the theoretical results 
of the decay to the ordinary p-shell states of n Be. With the help of the log ft values for 
the transitions from 11 Li, we can estimate the dissociation of a clusters in the daughter 
n Be nucleus. The observed strong GT transitions from n Li to n Be at 8 MeV comes from 
the breaking of the one of the 2 a clusters in the excited state of the daughter nucleus 
"Be. 
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